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Abstract 



o 

(N 

^0 , We explore the relation between fast waves, damping and imposed 

noise for different scalings by considering the singularly perturbed 
stochastic nonlinear wave equations 

, vutt + ut = Au + fiu) + iy'^W 

' on a bounded spatial domain. An asymptotic approximation to the 

r~| ! stochastic wave equation is constructed by a special transformation and 

splitting of vut ■ This splitting gives a clear description of the structure 
of u. The approximating model, for small > , is a stochastic non- 
linear heat equation for exponent < a < 1 , and is a deterministic 
nonlinear wave equation for exponent a > 1 . 



> 

^ ■ Keywords Singular perturbation, stochastic wave equations, asymptotic 

■ approximation. 

o 

cn 

0\ ■ 1 Introduction 

O . 

Our stochastic model is motivated by some material continuum in some 
domain D C M", 1 < n < 3 . The continuum is made of 'particles' with 
'displacement' field u{t,x) and 'velocity' field v{t,x). The motion of the 
^ ' particles in the continuum in a stochastic force field aW, motivated by 

■ Newton's law, is assumed to be described by the following stochastic partial 
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differential equations ^ 

v'it^x), (1) 
-kv^it, x) + Au^it, x) + fiu^it, x)) + aW{t, x), (2) 
uq , v^id.x) = ui , (3) 

for times t > , and locations x €z D , with zero Dirichlet boundary con- 
dition on dD. Here small u is the 'density' of the particles: we explore 
the singular limit as — t- and so label the dependent fields with super- 
script ly. Damping is proportional to the velocity v'^ with constant k. The 
Laplacian governs near neighbour, particle-particle, quasi-elastic interaction 
forces, and a nonlinear reaction is characterised by f{u). The externally 
imposed stochastic force field is (jW{t,x) where W{t,x) is an L2(L)) valued 
Wiener process defined on some complete probability space P), and 

is assumed to be of strength that scales according to a = i^". For exponent 
a = , the approximation of displacements n'^ , as — )• , is called the 
infinite dimensional Smolukowski-Kramers approximation which has been 
proved valid in the limit by estimating the remainder term [H [2] . Our re- 
cent work [5j applied an averaging method to approximate the displacement 
field u'^ for the case a = with exponent < a < 1/2 . Both of these 
methods are significantly complicated due to the coupling of displacement u 
and velocity v in the remainder term. 

Here we apply a relatively simple method to derive suitable approxima- 
tions for equations ([I])-© with a = i^", for exponent a S [0,1) U (l,oo). 
We apply the following useful splitting of the velocity u^, 

<{t) = U^At) + + i^'^-'^'vm, (4) 

to avoid directly estimating the remainder term. The three parts of the 
above splitting are the initial value part, the mean value part, and the 
diffusion part, respectively. This decomposition gives a clear structure for 
the displacement u'^ with 

^ Jo Jo Jo 

The parts and satisfy linear equations, and section [3] establishes 

- [ v'({s) ds = 0{u) and i/""^/^ / vlds = C(z/") as z/ ^ . 
^ Jo Jo 

The mean part (t) of the velocity is 0{1) as — )• for t G [0, T] with any 
fixed time T. Then for small section [3] determines which term is a high 
order term and gives an asymptotic approximation of the displacement vy . 

Here one interesting case is when the exponent a = 1 . In this case there 
are two terms with the same order 0{i') as — )• . Then if we keep all the 



vv^{t,x) = 
u^(0,3;) = 



2 



terms, the approximation to the displacement u'^ is just itself which is 
no modelling simplification. This case will be discussed further research. 

Because of its motivation by physical continuum problems of wave mo- 
tion in some random media [3], the system ([U-(l3|) is called a stochastic 
wave equation. For small ly and the particular case of o" = z^^/^ , Lv and 
Wang [6l [8] studied the limit behaviour as v ^ 0: in this case the random 
dynamics of ([I])-© was proved to be described by that of the nonlinear 
heat equation 

ut(t,x) = Au{t,x) + f{u{t,x)), u{0,x)=uo. (5) 

This paper extends this earlier research by approximating the behaviour 
of the solution on finite time interval [0,T], T > 0, for the more gen- 
eral case of fj = z^" with any a £ [0,l)U(l,oo). The tightness in the 
space C(0,T;L2(D)), compact in sense of probability, has been proved 
in previous work [21 |8]. Consequently, here we just need to approximate 
the displacement u'^ in a weak sense; that is, we consider the approxima- 
tion of the inner product {u'^,(p) in the space C(0,T) for testing function 
if G C^{D X [0,r]) with ip vanishing on the boundary dD. 

Section [2] first gives some preliminaries and the main result. Theorem HI 
Then section [3] details the proof. 



2 Preliminary 

Let D C M", 1 < n < 3 , be a regular domain with boundary dD. Denote 
by L?'{D) the Lebesgue space of square integrable real valued functions on D, 
which is a Hilbert space with inner product 

{u,v) = / u{x)v{x)dx, u,v£L'^{D). 
Jd 

Write the norm on L'^{D) by ||u||o = {u,u)^^^. Define the following abstract 
operator 

Au = -Au , u G Dom(^) = {n G L^{D) : Au G L'^{D) , u\dD = 0}. 

Denoted by {A^}, assume the eigenvalues of operator A satisfy < Ai < 
A2 < A3 < • • • , and Afc — )• 00 as /c — )• 00 . For any s > , denote by H^{D) the 
usual Sobolev space W'^'^iD) and by H'q{D) the closure of C^{D) in H'{D). 
In the space Hq{D) we use the equivalent norm 

\\u\\s = \\A'/^u\\o , u€H^{D). 

We also denote the dual space of Hq by H~^. Here specify that the noise 
magnitude scales as cr = z^", < < 1 , for exponent q > , in equation ([2]); 
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that is, we consider the following stochastic equations 

< = v\ n-(0) = no, (6) 
v'i = ^i-v" - Au'' + f{u'')] + u''~^W , <(0) = 'ui. (7) 

Hereafter we non-dimensionalise the time scale with the drag rate so that, 
in effect, the drag coefficient is one. We assume {W{t,x)}teM. is an L^{D)- 
valued, two sided, Wiener process, defined on a complete probability space 
J^, {J-t}t>o, P) with covariance operator Q such that 

Qek = hkCk , /c = 1,2, . . . , 

where {e^} is a complete orthonormal system in H and {h^} is a bounded 
sequence of non-negative real numbers. Then the noise process W{t,x) has 
the spectral expansion 

oo 

W{t,x) = ^ VhekWk{t), 

k=l 

where Wk are real, mutually independent, standard scalar Brownian mo- 
tions [3- Further, we assume boundedness of the sums 

oo oo 

trQ = bk < oo and Xkbk < oo . (8) 

fc=l k=l 

Assumption 1. For the nonlinearity f we assume 

1. 1/(^)1 <Ci(l + |5|3), |/'(S)|<C72(1 + |S|2); 

2. Fis) < -C3i\s\^ - 1), sfis) < -C^iFis) - 1); 

for some positive constants Ci, i = 1,2, 3, 4, and where F{s) = Jq f{r) dr . 

One simple example satisfying these assumptions is f{u) = u — . 

Then we have the following theorem. 
Theorem 2. Assume that the boundedness (0) and Assumption{l\hold. For 

any {uo,ui) € Hq{D)x L^[D) , there is a unique solution {u^ ,v^) to 
with 

u" e L\n,C{0,T;H}i{D))) and v" e L'^{n,C{0,T; L'^{D))), (9) 

for any T > . Moreover, for any T > there is a positive constant Ct 
which is independent of v such that the expectation 

E sup \\u''it)\\i<CTi\\uo\\l + \\ui\\l), 

0<t<T 

and {u^}o<i/<i is tight in the space C(0, T; L^(L')) . 
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Proof. To prove the existence of the solution we define 



A 



idi2 







and W{t) 





-1/2 



W{t) 



Let = {u'^,v'^), so equation ([S])-© can be rewritten in the following 
abstract stochastic evolutionary form 



6 = ^$ + F($) + >V, $(0) = (uo,ui). 



(10) 



Notice that operator A generates a strong continuous semigroup and the 
nonlinearity F is locally Lipschitz continuous, then by a standard method 
for stochastic evolutionary equations [7J we have the first part of the theorem. 

For < a < 1/2 the energy estimate for {u'^ , v^) and tightness result 
can be obtained via a similar argument to that of Cerrai and Freidlin [2], 
and for a > 1/2 the energy estimate and tightness were obtained by Lv and 
Wang |6j. The proof is complete. □ 

In the following approach we need the following lemma on weak conver- 
gence of a sequence of functions due to Lions 

Lemma 3. For any given functions hy and h G -L^([0,T] x D) {1 < p < 
oo), if ||/i!y||Lp([o,T]xD) ^ C for some positive constant C, and ^ h on 
[0, T]xD almost everywhere asv^Q, then h^ ^ h weakly in L^([0, T]xD). 

Now we give the main theorem on the approximation of the displace- 
ment vy in our stochastic wave equation. 

Theorem 4. Assume that the houndedness ^B^ and Assumption{I\hold, and 
(no, ui) £ Hq{D) X L?'{D) . If exponent < a < 1 , for any T > , and for 
small > , then with probability one 

U-'^Wu" - u''\\c(o^T;L\D)) ^ , aS U ^ , 

with the approximation solving the stochastic nonlinear heat equation 

< = Ar + /(r) + n'^(0) = no. (11) 

Conversely, if a > 1 , 

W'^Wu" - u''\\c{0,T;L^{D)) ^ , aS U ^ , 

with the approximation solving the deterministic nonlinear wave equation 



vutt + < = An'^ + /(n'^) , n"(0) = no , n^O) = ni . 



:i2) 
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3 Approximation: proof of Theorem [4] 

By Theorem [21 {u'^(t)}o<i/<i is tight in the space C(0, T; so we 

approximate the displacement in a weak sense: we approximate {u'^,ip) 
for any ip G C'^{[0, T] x D) with ip\gD = . 

In order to avoid the couphng between the displacement u'^ and the 
velocity v"^, we scale the velocity field as 



i^v". (13) 



Then 



< = 


V 




1 


= 


?} 

V 



u''(0) = no 



Further, we make the decomposition 

r = v\^vv'i^v''^^l'^vl, (14) 
where -u^^ = v\ ^ v\{^) = uu\^ (15) 

v^t = -hv2-^u''-f{un], vm = o, (16) 

Then 

= 1^5^ + t;^ + i/"-i/2^^ ^ ii^(0) = no. (18) 

The decomposition of v'^ makes the problem easier. The two sdes ()15p and ()17p 
for the two components ujj' and are just linear sdes whose properties are 
well known. The properties of V2 can be derived straightforwardly from the 
PDE ()16p by the estimates in Theorem [2j We state the following results. 
Lemma 5. Assume that the houndedness (0) and Assumption[J\ hold. Let 
ui G L'^{D), then for any ip G C^([0,T] x D) with ip\dD = 0, 

1 /■* 

- {v'({s),<fis))ds-^0, 0<t<T, (19) 
^ Jo 

and u-^''^ f {vl{s),ip{s))ds ^ [ {p{s),dW{s)), 0<t<T,{20) 
Jo Jo 

in L'^{fl), as u ^ . 

Proof. The proof is direct. First, 

vUt) = i^uie-'/". 
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Then for 99 G C'^{[0,T] x D) 



t/u 

(ui, Lp{vT))e~'^ (ir — )• , as — )• , 

uniformly on [0, T] , which yields the first convergence. Second, for any 
e C2([0,T] X D) with ^p\9o = by equation ([TT]), 

ly-'/' [\vUs),v{s))ds = -V^ f\vlt{s),^{s))ds+ [\ip{s),dW{s)) 
Jo Jo Jo 

Jo 

{^is),dWis)). (21) 



So it remains to show that vi^{t) is uniformly bounded in the space L'^{Q, L'^{D)). 
By equation p!7|) . applying Ito formula to H^g ||q gives 



Then by the Gronwall lemma 

E||v^(t)||o < trQ, t>0. 

The proof is complete. □ 

Lemma 6. Assume the conditions in Theorem holds, then there is a pa- 
rameter u and independent positive constant Ct such that 

^\v':^{t)\\-l < Ct , 0<t<T. 
Proof. For any ip G Hq{D) , from equation ()16p 

Then 

{vm,^) = -e'""" f e''" [-(Vn'^(s), VV^) + V)] ds . 



By the estimates in Theorem [2] and the embedding Hq{D) C L^{D) for 
1 < n < 3, 

E|(^2^(t),V)| <Ct||V||i, 0<t<T. 
The proof is complete. 

□ 
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Prom the above lemma we prove the main Theorem SI First, for any 
K > , by the tightness of displacement u'^ in the space C{0,T; L'^{D)), 
there is a compact set 5^ C C{0,T; L'^{D)) such that 

Flu" eB^}>l- Hi/2 . (22) 

By the Markov inequality and the estimate in Lemma [6l for any k > there 
is a positive constant such that 



"sv^^ll-i <Crll} > l-'^/2. (23) 
Then for any k > , define a probability space (fi^, J^k, IPk) 

n^ = {uj en-, events ([22]) and hold}, = {F n Q^, : F e T}, 
and for any F G J^^ 



P(fi«) 

In the following we restrict our problem to the above new probability space. 
For any u G 0,^ , the convergence still holds. The convergence ([^U]) is 
in the L^(Q) sense which yields the convergence for P almost all cj G , 
then we also have the convergence (|20p for almost all a; G $7^ . So we can 
assume that for all w G fi^ , the convergence ([20|) holds. 
Furthermore, we establish the limit 

f{u^)^f{u) weakly in L2(0,r;L2(D)) (24) 

for any ^ u in C(0, T; L2(i:))) . By the embedding H^{D) C L'^iD) we 
have ||/('u"')||L2{o,T;L2i:)) < Ct for some positive constant Ct , and by the 
strong convergence of n" — )• u , f{vP'{t, x)) — )• f{u{t, x)) on [0, T] x almost 
everywhere. Then Lemma [3] gives the limit. 

Next we give an asymptotic approximation to the displacement u'^ . For 
this we consider {u" {t) , ip{t)) with G C^{[0,T] x D) and lpIqd = 0. From 
equation (fT8]l 



K(t),<^(t))-(^zo,(^(0))- [ {u''{s),ipt{s))ds 

Jo 

= - [\v'as),ipis))ds+ f\v-^{s),^{s))ds + v^-^l^ f\vl{s),^{s))ds. 
^ Jo Jo Jo 

From equation 

\^;2'(s),^(s))^^5 = / {u''{s),Aip{s))ds+ / (/(n'^(s)),^(s))ds 
JO Jo 

- v{v2{t), ip{t)) + u I (7)2(5), ipt{s)) ds . 
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By the definition of i}^, and Lemmas [5] and [6l 



1 



V 

and 









Further, by 

l\vl^s),^{s))ds = u'^ f\^{s),dW{s))+0{u''+^/^). 
Jo Jo 

Then 

K(t),V^(t))-(no,¥'(0))- l\u%s),ipt{s))ds- f\u^{s),^^{s))ds 

Jo Jo 

{fiu^{s)),^is))ds 

= f\yi^(^s),^{s))ds + - [\v'[{s),^{s))ds 

Jo ^ Jo 

- iy{v2{t),ip{t)) + / {v2is),iptis)) ds 
Jo 

= [ {^{s),dW{s)) + 0{u"+^/^) + 0{u). (25) 
Jo 

Now for exponent < a < 1 , noticing the convergence (|24p and neglecting 
the o(i^") terms we have the following equation: 



(n"(t),(^(t))-(no,(/^(0))- / {u''{s),Ms))ds- / {u'' (s) , A^{s)) ds 

Jo Jo 

f ft 

\f{u%s)),^{s))ds = y'^ / {^{s),dW{s)). 

Jo 

Then we deduce the following approximation equation holds: 

v!i = Au" + fiu") + v^W , u^iO) = uo . (26) 

For exponent a > 1 we need to show the rate of decay in u of Vi and 
as 1/ — > . First, as {u'^}u C B,^ which is compact in C{0,T; L'^{D)) , for 
any sequence, there are a subsequence, say u'^" , with — )• as n — c« , 
and u £ C{0,T;L^{D)) with n(0) = uq such that 

u'^" — )• n , n — )• oo . 

First we assume 

Au + /(n) / in for all t G [0, T]. 
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Then by equation ([T6]), for any ip £ C'^i[0,T] x D) with ip\QD = 



Jo 

Jo Jo 

By the estimates in Lemma [6l and by the convergence ()24p . 

t f-t f-t 

{v!^"{s),ip{s))ds^ / {u{s),Aip{s))ds+ / {fiu{s)),ip{s))ds^O 
Jo Jo 

as n — )• oo . This shows that v^^it) does not converge to zero as — )• for 
<t<T. Then 

iy{v!^{t),ip{t)) =0{iy) and u [\v!^{s),^t{s)) ds = 0{iy) (27) 

JO 

instead of o(z^) as — )• . 

Further from equation (fT5]l , for any 99 G C^([0,T] x D) with (^I^d = , 

-KW, ^(i)) - {ui,m) - - fKis), Ms)) ds = -^ [\vUs), M)) ds . 
V V Jo ^ Jo 



Then 



1 



t 

(txi,(^(0)), i/^O. (28) 







Then for exponent a > 1 , in the asymptotic expansion ()25p for small 
neglecting the o(y) term consisting of U3, and by the transformation ()13[ 
and decomposition (JT 



(n'^(0,V5(t)) -(no,(/^(0))- [\u''{s),Ms))ds- ['{u^s), Aip{s)) ds 

Jo Jo 



V 



{fiu''{s)),^is))ds 



-/ (viis), ^{s)) ds - i/{v2{t), (fit)) + u I {v2{s),(ptis))ds 
Jo Jo 

i^{v''{t),M)) + u [ {v''{s),Ms))ds + i^{ui,ip{0)) 
Jo 





Then noticing that velocity v'^ = , and neglecting the 0(1/""^^/^) terms, 
we have the following approximation equation 

i/< + Ut = + /(r) . (29) 
10 



Second if for t G [0, T] 



Au + f{u)=0, u{0) = uo 

in H~^. Then it is a stationary solution of (|29p . 

The above approximation is in the sense of almost surely. Then 
by the arbitrary choice of k, and the well-posedness of (f26|) and ([29]) . this 
establishes the approximation with P probability one. 

This completes our proof of the main Theorem H] on the approximations 
of the stochastic wave equation ([I])-® for different scaling of the noise 
process. 

Acknowledgements This research was supported by the NSFC grant 
No. 10901083 and by the Austrahan Research Council grants DP0774311 
and DP0988738. 

References 

[1] S. Cerrai & M. Freidlin, On the Smoluchowski-Kramers approximation 
for a system with an infinite number of degrees of freedom, Prob. Th. 
and Relat. Fields 135 (2006), 363-394. 

[2] S. Cerrai & M. Freidlin, Smoluchowski-Kramers approximation for a 
general class of SPDEs, J. Evol. Equa. 6 (2006), 657-689. 

[3] P. L. Chow, W. Kohler, & G. Papanicolaou, Multiple Scattering and 
Waves in Random Media, North-Holland, Amsterdam, 1981. 

[4] J. L. Lions, Quelques methodes de resolution des problemes non 
lineaires, Dunod, Paris, 1969. 

[5] Yan Lv & A. J. Roberts, Averaging approximation to singu- 
larly perturbed nonlinear stochastic wave equations, preprint, 2011. 
|http : //arxiv^ org7abs/l 107 . 4184 

[6] Y. Lv & W. Wang, Limit dynamics for stochastic wave equations, J. 
Diff. Equa. 244 (2008), 1-23. 

[7] G. Da Prato & J. Zabczyk, Stochastic Equations in Infinite Dimensions, 
Cambridge University Press, 1992. 

[8] W. Wang and Y. Lv, Limit behavior of nonlinear stochastic wave equa- 
tions with singular perturbation. Disc, and Cont. Dyna. Syst. B, 13(1) 
(2010) 175-193. 



11 



